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Abstract 

We calculate the full electroweak one-loop corrections to the decay of the CP- 
odd Higgs boson A into scalar quarks in the minimal supersymmetric extension 
of the Standard Model (MSSM). Due to the complex structure of the electroweak 
sector a proper renormalization of many parameters in the on-shell renormalization 
scheme is necessary. For the decay into sbottom quarks, especially for large tan/3, 
the corrections can be very large in the on-shell renormalization scheme, which 
makes the perturbation series unreliable. We solve this problem by an appropriate 
definition of the tree-level coupling in terms of running quark masses and running 
trilinear couplings A q . We also discuss the decay of heavy scalar quarks into light 
scalar quarks and ^4°. We find that the corrections are significant and therefore 
cannot be neglected. 



1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) pQ requires five physical Higgs 
bosons: two neutral CP-even (h° and H°), one heavy neutral CP-odd (A ), and two 
charged ones (H^) [3 El- The existence of a CP-odd neutral Higgs boson would provide 
a conclusive evidence of physics beyond the SM. Searching for Higgs bosons is one of the 
main goals of present and future collider experiments at TEVATRON, LHC or an e + e~ 
Linear Collider. 

In this paper, we consider the decay of the CP-odd Higgs boson A into two scalar 
quarks, A — > q\q 2 . The decays into squarks can be the dominant decay modes of Higgs 
bosons in a large parameter region if the squarks are relatively light jHE]. In particular, 
the third generation squarks t« and fej can be much lighter than the other squarks due 
to their large Yukawa couplings and their large left-right mixing. We will calculate the 
full electroweak corrections in the on-shell scheme and will implement the SUSY-QCD 
corrections which were calculated previously j^j. The challenge of this calculation is the 
necessity to renormalize almost all parameters in the electroweak sector in only one single 
process. Due to the numerous electroweak interacting particles and the complex coupling 
structure we have to compute a large number of graphs. In general, the Higgs-squark- 
squark couplings consist of F— and D-terms and SUSY breaking terms, all depending on 
the squark mixing angle 9 q . As a first step we consider the case A — > qi q 2 where only 
F-terms and SUSY breaking terms enter in the coupling. Since A only couples to qi^qn 
and due to the CP nature of A , A — > q~i q i vanishes (with real parameters also beyond 
the tree-level!). Despite the complexity, we have performed the calculation in an analytic 
way. The explicit formulae will be given elsewhere. We will, however, show the most 
important results of the numerical analysis. Furthermore, the crossed channel q 2 -> qtA 
is studied. 

In case of the decay into sbottom quarks the decay width can receive large corrections 
which makes the perturbation expansion unreliable, especially for large tan/3. In some 
cases the width can even become negative using the on-shell renormalization scheme. 
We will show that this problem can be fixed by an appropriate choice of the tree-level 
coupling in terms of DR running quark masses and running A q . 

2 Tree— level result 

The squark mixing is described by the squark mass matrix in the left-right basis (q L , q R ), 
and in the mass basis (q~i, q 2 ), q = t or b, 




(1) 



where Rf a is a 2 x 2 rotation matrix with rotation angle 9 q , which relates the mass 
eigenstates q^ i — 1, 2, (m qi < m q2 ) to the gauge eigenstates q a , a = L,R, by q^ = Rl a q a 
and 



ml = Mq + (Iq L — e q sm 2 9 w ) cos 2(3 m\ + m 2 q , (2) 
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m qR = M0 m + e <? sin ®w cos 2(3 m z +m q , (3) 



A q - ii (tan /3)~ 2I * L . (4) 



Mq, Mjj, and are soft SUSY breaking masses, A q is the trilinear scalar coupling pa- 
rameter, n the higgsino mass parameter, tan f3 — ^ is the ratio of the vacuum expectation 
values of the two neutral Higgs doublet states 0] , -ff X denotes the third component of 
the weak isospin of the quark q, e q the electric charge in terms of the elementary charge 
e , and 6 W is the Weinberg angle. 

The mass eigenvalues and the mixing angle in terms of primary parameters are 



m L = o \ m \, + m L =f \l ( m L ~ m L) 2 + Aa l m 1 ) ( 5 ) 
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cos^~ = 9 g (0<%<tt), (6) 

and the trilinear breaking parameter A q can be written as 

m q A q = - (m| - m~ 2 ) sin 26» 9 - + m q /j (tan/3)~~ 2/ ' L . (7) 
At tree-level the decay width of A — * q~i q 2 is given by 



r tree M C> 3 K Wo » m L m | ) | 

1 [A qiq 2 ) — — 3 |<^i23 



A 

with y, z) = J(x — y — z) 2 — Ayz and the A°-q*-qj coupling Gfj 3 given in [§]. 



3 Full Electroweak Corrections 

The one-loop corrected (renormalized) amplitude G123 11 can be expressed as 

^123 — ^123 + ^^123 — ^123 + 0Lr l23 + 0ljr 123 + °^123 ; l y J 

where ^Crf^ are the vertex corrections (Fig. and SGl^ the wave-function corrections 
(Fig.|2J). Note that in addition to the one-particle irreducible vertex graphs also one-loop 
induced reducible graphs with A°-Z° mixing have to be included. All parameters in the 
tree-level coupling Gi 23 have to be renormalized due to the shift from the bare to the 
on-shell values. These corrections are denoted by <?G©. The full one-loop corrected 
decay width is then given by 

r(A° - ~ qi q 2 ) = *« m *><><) [| G f 23 | 2 + 2Re ( G f 23 . AGfa3 )] . (10) 

lbirm A Q 1 v /J 

Since there are diagrams with photon exchange we also have to consider corrections due to 
real photon emission to cancel the infrared divergences (Fig. Therefore the corrected 
(UV- and IR-convergent) decay width is 

pcorr^O _^ ~ = T{A ^ - + p ( ^0 ^ ~ (n) 
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Throughout the paper we use the SUSY invariant dimensional reduction (DR) as regu- 
larization scheme. For convenience we perform the calculation in the 't Hooft-Feynman 
gauge, £ = 1. 



3.1 Vertex and wave— function corrections 

The relations between the unrenormalized (bare) and renormalized (physical) fields and 
couplings are 

C n = C TC11 + 5C, 



Co = W W° • (Gf 23 y = G\ 23 + 5G 

£ ren = G\ 23 A ~q\q 2 , (A°) ° = y/l + 5Zg Hi , 



123 5 







^A%q 2 , (qt) =^l + 5Ziq* 

, o 



(12) 



92 =V 1 + ^2 J ^ 



with the notation = {h°, H°, A , G }, i,j = 1, 2, and k = 3,4. Thus the renormalized 
Lagrangian is given by (up to the first order) 

£ ren = (CL + SGS + \ {SZlGf 23 + SZ] 2 G\ j3 + 6Z&GQ + 5Gfg) A°q{q 2 . 

(13) 

The explicit form of the vertex corrections 5Gf 23 will be given elsewhere. Due to the 
antisymmetry of the tree-level coupling, Gfj 3 = —Gj i3 , the total wave-function correction 
reads 

SGS = l(5Z! 1 + 5Z! 2 + 5Zg)Gt 23 + l-5ZgGl 24 . (14) 



2 \ ii ' zz 1 j ±zo ' 2 

For the wave-function renormalization constants we use the conventional on-shell renor- 
malization conditions [2j which lead to 

5Zl = -»nj(mj,) , 2 ff/ 2 N 

6Zg = — -^Uf 3 (m%), 15 

*Z* = -»nS(mio), m 2 G o-mio 

with the diagonal parts of the Higgs and squark self-energies Ylu(k 2 ). 
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Figure 1: Vertex and photon emission diagrams relevant to the calculation of the virtual 
electroweak corrections to the decay width A — > <?i <? 2 - 
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Figure 2: Wave-function diagrams relevant to the calculation of the virtual electroweak 
corrections to the decay width A — > q± q 2 . H k denotes neutral as well as charged Higgs 
bosons. 

The off-diagonal Higgs wave-function corrections can be combined with the contribution 
to which come from A°-Z° mixing. First we show that the sum of the parts 

coming from the propagators of Z° and G° outside the loops is independent of the gauge 
parameter £ = 
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Figure 3: A°-Z° contribution and A°-G° wave-function correction 
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The amplitudes of the two graphs of Fig. El in a general gauge are 



M z = (-ip^U AZ (p 2 



p 2 —m 2 



-<^+(W) 



p 2 —£m 



Wz 42)(ki+k 2 ) l/ 



Contracting the Lorentz indices in Ai z , 

P^-g^ + il-O^P^Aih+hY 
V p l —t > m l z ) 



- 1 



p 2 — £m 



[ m qi - m h) 



and eliminating Hag in favor of Haz by using the Slavnov-Taylor identity jH] 

p 2 U A z(p 2 )+im z U A G{p 2 ) = 0, 
we find the sum M z + M G 



M 



Z+G 



p 2 — m 2 



n A z(p 2 )gzzl 2 (m] a ^ 



91 92 



p 2 — £m\ 



+ 



Finally we use the identity 



to obtain the result 



p 2 U AZ (p 2 



p 2 — £m| m 2 
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im z G q ij4 



tM z+G (p 2 - m^o) 



r 123 



ilU^m 2 ^) G'i24 



m z (p 2 — m|) (p 2 — £,rn 2 ) 



x 
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Ti A z(rn 2 AQ ) G\ 2A . 



(16) 
(17) 



(19) 



(20) 



f21] 



(22) 



The gauge dependence of the propagators of the Z° and G° in Figgis completely removed. 
However, there still remain gauge dependences from vector particles and Goldstone bosons 
in the loops of Haz which cancel against their counter parts in the vertex, wave-function 
and counter term corrections. 



3.2 Counter terms 

Since all parameters in the tree-level coupling G\ 2 s have to be renormalized, we get 



5G 



123 



Sh 
~h, 



cos 



^ 23 + ^ M K|— 



(23) 



for j d( ^ n }~tyP e squarks. The Yukawa coupling counter term can be decomposed into 
corrections to the electroweak coupling g, the masses of the quark q and the gauge boson 
W and the mixing angle /3, 

5h q 5g ^ 5m q 5m w ^ j - cos 2 (3 j 8 tan (3 
h q g m q m w \ sin (3 J tan (3 



For the trilinear coupling we get with eq. 

5A q S(m q A g ) 5m q 



A q TYlqAq TTlq 



(25) 



8(m q A q ) = 5 [rriqfi j Ij + ^(Sm^-Srn^) sin 20 q 



+ (ml-ml)cos29qS9q. (26) 

In the on-shell scheme the renormalization condition for the electroweak gauge boson 
sector reads [H] 

= — + 1 ( Sm w _ Sm z \ ^ 27 N 
g e tan 2 9 W \ m w m z J 

with m w and m z fixed as well as the quark and squark masses as the physical (pole) 
masses. 

Renormalization of the electric charge e 

Since we use as input parameter for a the MS value at the Z-pole, a = Oi{m z )\y^ = 
e 2 /(47r), we get the counter term [TU] 



5e 1 e 



e (4tt) 2 6 



4E^e 2 / (A + log^)+Y:E^4( A + lo g^) 

■ / V X f J f m=l V m fJ 

-4 E (A + log - E (A + log X-) - 2 (A + log 

j^V m 2 + J £\ m 2 + J V m 2 w J 



(28) 

with Xf = m z \/ rrif < m z and x t = m t . Nq is the colour factor, N[, = 1, 3 for (s)leptons 
and (s)quarks, respectively. A denotes the UV divergence factor, A = 2/e — 7 + log47r. 

Renormalization of tan (3 

For tan/3 we use the condition ImII^o z o(m 2 1 ) = which gives the counter term 

<5tan/3 1 T „ / 9 . . , 

— = — ImrW m 2 o . (29) 

tanp m z sm zp 
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Renormalization of /x 

The higgsino mass parameter /x is renormalized in the chargino sector ^21 02! where it 
enters in the 22-element of the chargino mass matrix X, 



Renormalization of 0q 



The counter term of the squark mixing angle, 56q, is fixed such that it cancels the anti- 
hermitian part of the squark wave-function corrections [Tlj I15j. 



59, 



•i 



\ (5Zf 2 - 5Zl) = 2(m /_ m ^ Re (nf 2 (m|) + 14^)) . (31) 



3.3 Infrared divergences 

The infrared divergences in eq. (fTU|) are cancelled by the inclusion of real photon emission, 
see the last two Feynman diagrams of Fig. [T] The decay width of A°(p) — > qi(ki) + q 2 (k2) + 
7(^3) can be written as 

r(A°-M 27 ) = 3( ^ g) 3 |Gf23 ' 2 [(m%-ml-ml) h,-m\ x h x -m\h 2 -h-h 

(32) 

with the phase-space integrals I n and I mn defined as [HI] 

1 f^hd^d 3 ^^ 1 

h —--zJ mm'm {p 1 2 3j (2^ + a 2 ) . . . (2^ + a*) ■ {66) 

The corrected (UV- and IR-convergent) decay width is then given by (see eq. ([11))) 

pcorr^O ^ = r(A ^ + ^0 ^ - ^ ^ _ ^ 



4 Improvement of One— loop Corrections 

In the on-shell renormalization scheme, in case of the decay into sbottom quarks, es- 
pecially for large tan f3, the decay width can receive large corrections which makes the 
perturbation expansion unreliable. In some cases the corrected width can even become 
negative. It has been pointed out [HI HHj that the source of these large corrections are 
mainly the counter terms for m& and the trilinear coupling A&. We show that this problem 
can be fixed by absorbing these large counter terms into the y4°-squark-squark tree-level 
coupling and expanding the perturbation series around the new tree-level. The technical 
details will be given in a forthcoming paper. 
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Correction to m.5 

If the Yukawa coupling h b is given at tree-level in terms of the pole mass mj, the one- 
loop corrections to the counter term Sm b become very large due to gluon and gluino 
exchange contributions. We absorb these large counter terms and also the ones due 
to loops with electroweak interacting particles into the Higgs-squark-squark tree-level 
coupling by using the DR running mass rh b (Q = rriA)- The large counter term due to the 
gluon loop is absorbed by using SM 2-loop renormalization group equations [THl EE 120] • 
Thus we obtain the SM running bottom m b (Q)sM- For large tan j3 the counter term to m b 
can be very large due to the gluino-mediated graph [T71 |2U 122] • Here we absorb the gluino 
contribution as well as the sizeable contributions from neutralino and chargino loops and 
the remaining electroweak self-energies into the Higgs-squark-squark tree-level coupling. 
In such a way we obtain the full DR running bottom quark mass 

m b (Q) M ssM = m b (Q) SM + 5m b (Q) . (35) 

Correction to Ab 

The second source of a very large correction (in the on-shell scheme) is the counter term 
for the trilinear coupling A b , eqs. (I2H| EEJ), especially the contribution of the left-right 
mixing elements of the squark mass matrix, m 2 LR = (m? — m? 2 ) sin^cos^q. As in the 
case of the large correction to m b we use DR running A^m^o) in the Higgs-squark-squark 
tree-level coupling. Because of the fact that the counter term 5A b (for large tan/3) can 
become several orders of magnitude larger than the on-shell A b we use A^m^o) as input 
[T%] . In order to be consistent we have to perform an iteration procedure to get the correct 
running and on-shell masses, mixing angles and other parameters. 

5 Numerical analysis and conclusions 

In the following numerical examples, we assume Ma = Ma = yM^ 3 = j^M^ = Mi 3 = 
M Es = M Qi 2 = M Di2 = M Di2 = M Li 2 = M£ 12 for the first, second and third gen- 
eration soft SUSY breaking masses and A = A t = A& — A T , if not stated otherwise. 
For the standard model parameters we take m z = 91.1876 GeV, m w = 80.423 GeV, 

sin 2 ^ = 1 - m^/ml, a = 1/127.934, m t = 174.3 GeV, and m b = 4.7 GeV. W is fixed 

5 

by the gaugino unification relation M' = — tan 2 9 W M and the gluino mass is related to 
M by nig = (a s (m g )/a) sin 2 9 W M. 

Decays into stops: 

In Fig.0]we show the tree-level and the corrected width to A — > t\t2 for tan (3 = 1 and 
{Ma, A, M, fi} = {300, —500, 120, —260} GeV as a function of the mass of the decaying 
Higgs boson, m^o. As can be seen for larger values of m^o, the electroweak corrections 
can be of the same size as the SUSY-QCD corrections. 

In Fig. El the tree-level, the full electroweak and the full one-loop corrected (elec- 
troweak and SUSY-QCD) decay width of A -> tj 2 are given as a function of A t . The 
electroweak corrections do not strongly depend on the parameter A t and are almost con- 
stant about 8%. As input parameters we have chosen the values given above as well as 
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Figure 4: Tree-level (dotted line), full electroweak corrected (dashed line) and full one- 
loop (electroweak and SUSY-QCD) corrected (solid line) decay width of A —>■ t 1 t 2 - 
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Figure 5: ^-dependence of tree-level (dotted line), full electroweak corrected (dashed 
line) and full one-loop (electroweak and SUSY-QCD) corrected (solid line) decay width 
of A — > tit 2 . The gray area is excluded by experimental bounds. 



{A b>T ,m A o} = {-500,700} GeV. 

Fig. El shows the tree-level, the full electroweak and the full one-loop corrected (elec- 
troweak and SUSY-QCD) decay width of A — > t\t 2 as a function of tan (3 with the same 
parameter set as above and m A o = 900 GeV. Again, in a large region of the parameter 
space the electroweak corrections are comparable to the SUSY-QCD ones. 

Decays into sbottoms: 

Here we illustrate the numerical improvement of the full one-loop corrections to A —>■ &!& 2 
for large tan (3. 

In Fig. [7| we show two kinds of perturbation expansion for the input parameters 
{m A o, Mq, A u A b , At, M, fi} = {800, 300, 150, -700, -500, 120, 260} GeV: First we show 
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tan (5 

Figure 6: Tree-level (dotted line), full electroweak corrected decay width (dashed line) 
and full one-loop (electroweak and SUSY-QCD) corrected width (solid line) of A — > tit^ 
function of tan (5. 
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Figure 7: Two kinds of perturbation expansion: the dotted line corresponds to the on- 
shell tree-level width, the dashed and dash-dot-dotted line correspond to electroweak 
SUSY-QCD on-shell one-loop width, respectively. The dash-dotted line corresponds to 
the improved tree-level and the solid line to the (full) improved one-loop width. 



the on-shell tree-level width (dotted line). The dashed and dash-dot-dotted line corre- 
spond to the on-shell electroweak and full (electroweak plus SUSY-QCD) one-loop width, 
respectively. For both corrections one can clearly see the invalidity of the on-shell pertur- 
bation expansion, in particular the electroweak corrections lead to an improper negative 
decay width. The second way of perturbation expansion is given by the dash-dotted and 
the solid line which correspond to the improved tree-level and improved full one-loop 
decay width, respectively. The smallness of the relative correction in this case shows 
that the improved tree-level is already a good approximation for A — > h\\)2- The input 
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parameters are the same as in the first case but now with running A& = —700 GeV. 
Squarks decays: 

Fig. IS] displays the decay widths of the crossed channel t 2 — > t\A° as a function of 
A t . As can be seen, the electroweak corrections are as large as the SUSY-QCD ones in 
the considered region. The values of the input parameters are {tan/3, /i} = {35, —300} 
and {m A o, m- g , Mq, A b , A T } = {150, 1000, 300, -700, -700} GeV with the relations for the 
SUSY breaking masses given at the top of this section but with M^ 3 = 500 GeV in order 
to get a quite acceptable mass splitting in the stop sector. 



> 

o 

a 



o 




Figure 8: ^-dependence of the tree-level (dotted line), full electroweak corrected (dashed 
line) and full one-loop corrected (solid line) decay widths of t 2 — > tiA°. 
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Figure 9: Decay widths of b 2 — > biA° as a function of tan/3. The dotted and dash-dot- 
dotted line correspond to the on-shell tree-level and on-shell one-loop width, respectively. 
The dash-dotted line corresponds to the full improved tree-level and the solid line to the 
full improved one-loop width. 
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Fig. El again demonstrates the numerical improvement in the large tan/3 regime: The 
dotted and dash-dot-dotted lines correspond to the on-shell tree-level and on-shell one- 
loop decay widths of b 2 — > b\A°, whereas the dash-dotted and solid lines show the full 
improved tree-level and one-loop widths, respectively. The input parameters are the 
same as in Fig. H but with {Mq 3 ,A} = {500, -700} GeV. 

In conclusion, we have calculated the full electroweak one-loop corrections to the decay 
widths A — > qi q 2 and q 2 — > q\A° in the on-shell scheme. Moreover, we have included the 
SUSY-QCD corrections which were calculated in For the decay into sbottom quarks 
and large tan (3 an improvement of the on-shell perturbation expansion is necessary. This 
was done by an appropriate redefinition of the tree-level Higgs-squark-squark coupling. 
We find that the corrections are significant and in a wide range of the parameter space 
comparable to the SUSY-QCD corrections. 

Acknowledgements 

The authors acknowledge support from EU under the HPRN-CT-2000-00149 network 
programme and the "Fonds zur Forderung der wissenschaftlichen Forschung" of Austria, 
project No. P13139-PHY. 

References 

[1] H. E. Haber and G. L. Kane, Phys. Rep. 117 (1985) 75. 

[2] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, The Higgs Hunter's Guide, 
Addison- Wesley (1990). 

[3] J. F. Gunion, H. E. Haber, Nucl. Phys B 272 (1986) 1; B 402 (1993) 567 (E). 

[4] A. Bartl, K. Hidaka, Y. Kizukuri, T. Kon and W. Majerotto, Phys. Lett 315 (1993) 
360; 

A. Djouadi, J. Kalinowski, P. Ohmann, P.M. Zerwas, Z. Phys. C 74 (1997) 93. 

[5] A. Bartl, H.Eberl, K. Hidaka, T. Kon W. Majerotto and Y. Yamada, Phys. Lett 378 
(1996) 167 and references therein. 

[6] A. Bartl, H. Eberl, K. Hidaka, T. Kon, W. Majerotto and Y. Yamada, Phys. Lett. 
B 402 (1997) 303; 

A. Arhrib, A. Djouadi, W. Hollik, C. Jiinger, Phys. Rev. D 57 (1998) 5860. 

[7] M. Bohm, W. Hollik and H. Spiesberger, Fortschr. Phys. 34 (1986) 687. 

[8] A. Dabelstein, Z. Phys. C 67, 495 (1995); Nucl. Phys. B 456, 25 (1995). 

[9] A. Sirlin, Phys. Rev. D 22 (1980) 971; 

W.J. Marciano and A. Sirlin, Phys. Rev. D 22 (1980) 2695; 
A. Sirlin and W.J. Marciano, Nucl. Phys. B 189 (1981) 442. 

[10] H. Eberl, M. Kincel, W. Majerotto and Y. Yamada, Nucl. Phys. B 625 (2002) 372. 



14 



[11] P. H. Chankowski, S. Pokorski, J. Rosiek, Phys. Lett. B 274 (1992) 191; Nucl. Phys. 
B 423 (1994) 437; 497; 

A. Dabelstein, Z. Phys. C 67 (1995) 495; Nucl. Phys. B 456 (1995) 25. 

[12] H. Eberl, M. Kincel, W. Majerotto and Y. Yamada, Phys. Rev. D 64 (2001) 115013. 

[13] W. Oiler, H. Eberl, W. Majerotto and C. Weber, |hep-ph/0304"006l to be published 
in Eur. Phys. J. C. 

[14] J. Guasch and J. Sola, W. Hollik, Phys. Lett. B 437 (1998) 88. 

[15] H. Eberl, S. Kraml and W. Majerotto, JHEP 9905 (1999) 016. 

[16] A. Denner, Fortschr. Phys. 41 (1993) 307. 

[17] R. Hempfling, Phys. Rev. D 49 (1994) 6168; 

L.J. Hall, R. Rattazzi, and U. Sarid, Phys. Rev. D 50 (1994) 7048; 

M. Carena, M. Olechowski, S. Pokorski, and C.E.M. Wagner, Nucl. Phys. B 426 

(1994) 269; 

D.M. Pierce, J. A. Bagger, K. Matchev, and R. Zhang, Nucl. Phys. B 491 (1997) 3. 

[18] H. Eberl, K. Hidaka, S. Kraml, W. Majerotto and Y. Yamada, Phys. Rev. D 62 
(2000) 055006. 

[19] E. Braaten and J. P. Leveille, Phys. Rev. D 22 (1980) 715; 

M. Drees and K. Hikasa, Phys. Lett. B 240 (1990) 455; 262 (1991) 497 (E); 
A. Mendez and A. Pomarol, Phys. Lett. B 252 (1990) 461. 

[20] S.G. Gorishny, A.L. Kataev, S.A. Larin, and L.R. Surguladze, Mod. Phys. Lett. A5 
(1990) 2703; Phys. Rev. D 43 (1991) 1633; 

A. Djouadi, M. Spira, and P.M. Zerwas, Z. Phys. C 70 (1996) 427; 

A. Djouadi, J. Kalinowski, and M. Spira, Comput. Phys. Commun. 108 (1998) 56; 

M. Spira, Fortschr. Phys. 46 (1998) 203. 

[21] M. Carena, S. Mrenna, and C.E.M. Wagner, Phys. Rev. D 60 (1999) 075010; Phys. 
Rev. D 62 (2000) 055008. 

[22] P.H. Chankowski and S. Pokorski, |hep-ph/ 9707497, in Perspectives on Supersymme- 
try, ed. by G. L. Kane (World Scientific, 1997); 
K.S. Babu and C. Kolda, Phys. Lett. B 451 (1999) 77. 



15 



